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I  INTRODUCTION 


The  objective  of  this  project  was  to  research  NASTRAN's  effectiveness 
in  analyzing  nuclear  blast  overpressure  effects  on  panels  as  simulated  by 
shock  tube  tests.  Ultimately,  this  determines  NASTRAN's  effectiveness  in 
predicting  sure  safe  panel  response  to  nuclear  blast  overpressure  effects 
for  survivability /vulnerability  analysis.  Accomplishment  of  this  objective 
was  achieved  by  comparing  NASTRAN  data  to  experimental  shock  tube  test  data 
which  the  Defense  Nuclear  Agency  (DNA)  collected  with  the  Boeing  Military 
Airplane  Company  under  contract  DNA-OOI-76-C-0084  and  published  in  DNA 
report  DNA-4278F,  Volumes  1  through  4  (hereinafter  addressed  as  reference 
1).  Shock  tube  tests  are  an  accepted  method  for  simulating  the  effects  of 
nuclear  blast  overpressures.  Therefore,  comparing  NASTRAN  data  to  shock 
tube  test  data  is  an  effective  method  for  validating  NASTRAN  as  an 
overpressure  analysis  technique. 

NASTRAN  is  a  finite  element  structural  analysis  computer  code  that  is 
universally  accepted  in  the  structural  analysis  community.  The  version  of 
NASTkAN  used  in  the  analysis  for  this  report  is  COSMIC,  a  linear  analysis 
valid  only  for  predicting  panel  response  to  the  yield  point.  Experimental 
data  used  for  comparison  with  NASTRAN  came  from  shock  tube  tests  performed 
upon  seven  panel  configurations.  These  configurations  varied  in  thickness, 
edge  support  constraints,  magnitudes  of  subjected  overpressures,  geometry, 
and  materials  (see  Table  2  and  Figures  13  through  15).  Magnitude  of  first 
deflection  was  the  criterion  used  to  measure  NASTRAN'S  effectiveness. 

Stress  was  not  used  as  a  criterion  because  deflection  data  in  reference  one 


Is  of  more  consistent  quality.  Deflection  and  stress  exhibit  a  linear  rela¬ 
tionship  In  a  material's  elastic  range.  Therefore,  deflection  Is  a  valid 
measure  of  NASTRAN's  effectiveness  for  predicting  sure  safe  panel  response. 
This  Is  discussed  In  further  detail  In  Section  IV,  Discussion. 

Shock  tube  test  are  performed  by  generating  a  shock  wave  that  propaga¬ 
tes  down  a  tube  and  strikes  a  specimen.  Experimental  data  used  for  com¬ 
parison  In  this  project  was  performed  at  Sandia  Corporation's  THUNDERPIPE 
shock  tube  In  Albuquerque,  New  Mexico.  The  THUNDERPIPE  shock  tube  genera¬ 
tes  a  shock  wave  by  prlmacord  explosives.  Figure  1  is  taken  from  DNA  report 
DNA-4278F  and  Illustrates  the  dimensions  of  the  Thunderpipe  Shock  Tube.  The 
reader  should  realize  that  this  Is  a  relatively  large  test  facility. 
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II  SHOCK  TUBE  TEST  DATA 


Researching  NASTRAN's  effectiveness  for  reproducing  structural  respon¬ 
ses  observed  in  shock  tube  tests  required  development  of  the  data  interpre¬ 
tation  methodology  Introduced  in  this  section.  This  development  is  divided 
into  the  two  subsections:  Data  Interpretation  and  Error  Effects.  Large 
experimental  data  fluxuations  required  development  of  a  data  interpretation 
methodology.  This  methodology  provides  a  consistent  method  of  interpreting 
the  pressure  time  history  data  reported  in  reference  1.  The  interpreted 
data  is  Input  into  the  NASTRAN  model  built  to  simulate  the  tested  structure 
Error  Effects  is  a  study  of  the  effects  upon  NASTRAN  analysis  if  data  con¬ 
tains  an  Inherent  interpretational  error. 


II-l  DATA  INTERPRETATION 


Validating  any  numerical  analysis  technique  requires  accurate  and  con¬ 
sistent  methods  for  reading  experimental  data  used  as  analytical  input  data. 
The  following  is  a  methodology  developed  for  reading  experimental  data  of 
pressure  time  histories  produced  in  the  THUNDERPIPE  shock  tube.  Ideal 
overpressure  curves  for  reflected  pressure  time  histories  are  the  guidelines 
for  data  interpretations. 

The  reference  used  for  ideal  blast  waves  is:  The  Effects  of  Nuclear 
Weapons,  compiled  and  edited  by  Samuel  Glasstone  and  Philip  J.  Dolan,  3rd 
edition,  published  by  the  United  States  Department  of  Defense  and  the  Energy 
Research  and  Development  Administration  (hereinafter  addressed  as  reference  2). 
Ideal  curves  for  blast  waves  seldom  correlate  exactly  to  experimental  shock 
tube  data.  Thus,  it  is  emphasized  that  ideal  curves  are  used  only  as  guide¬ 
lines.  The  methodology  developed  pertains  to  ideal  curves  for  surface  blast 
waves  that  strike  normal  to  flat  and  curved  panels.  Section  II-l.A, 

Definitions  of  Terms,  will  enhance  the  reader's  understanding  of  the  metho¬ 
dology. 

Reflected  pressure  spikes  are  the  most  important  consideration  when 
interpreting  pressure  time  history  data.  Spike  peaks  and  widths  are  depen¬ 
dent  upon  post-reflected  peaks.  Therefore,  interpreting  experimental 
reflected  pressure  data  requires  that  post-reflected  curves  be  determined 
first,  followed  by  interpretation  of  spike  peaks  and  then  spike  peak  widths. 
Fitted  curves  will  vary  from  one  interpreter  to  another,  but  the  differences 
will  be  negligible  if  the  guidelines  for  this  methodology  are  followed. 
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II-1.A  DEFINITIONS  OF  TERMS 


Reflected  Pressure  (Pr)  -  The  Initial  pressure  experienced  by  surfaces  sub¬ 
jected  to  shock  waves  striking  at  non-parallel  paths  relative  to  the 
surface,  resulting  In  a  greats  -  pressure  experienced  by  the  surface 
than  is  present  at  the  shock  front. 

Reflected  Pressure  Time  History  -  A  numerical  account  (tabular  or  graphical) 
of  the  pressure  as  a  function  of  time  experienced  by  a  surface  subjected 
to  shock  waves  traveling  non-parallel  paths  relative  to  the  surface.  It 
Is  the  addition  of  incident  overpressure,  dynamic  pressure,  and 
reflected  pressure  effects. 

Stagnation  Time  (ts)  -  The  time  at  which  reflected  pressure  effects  subside, 
leaving  only  incident  overpressure  and  dynamic  pressure  effects.  It  is 
a  function  of  panel  geometry  and  shock  wave  velocity.  Ideal  ts  is 
calculated  from  reference  2. 

Stagnation  Pressure  (Ps)  -  The  post- reflected  peak  pressure  that  corresponds 
to  stagnation  time  (ts). 

Ideal  Pressure  Curve  -  Developed  in  reference  2.  Major  characteristics  are 
an  initial  reflected  pressure  effect  until  time  ts,  followed  by  a 
steady  and  more  gradual  decrease  of  the  post-reflected  pressure. 


SPIKE 

PEAK- 

WIDTM 


REFLECTED 
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ZONE 


II-l.B  POST-REFLECTED  CURVES 


Experimental  shock  wave  data  exhibits  extreme  fluxuations  in  the  post- 
reflected  pressure  zone  (see  Figure  3).  These  fluxuations  are  due  to  the 
combined  effects  of  multiple  detonations  during  ignition  of  the  primacord, 
close  proximity  of  the  test  specimen  to  the  explosive,  and  possible  experi¬ 
mental  data  noise.  A  realistic  approach  to  analysis  requires  this  data  be 
approximated  as  a  smooth  curve.  Time  steps  required  to  analyse  actual 
experimental  data  fluxuations  would  result  in  unnecessary  expenditure  of 
computer  time.  Selecting  the  post-reflected  curve  is  accomplished  by 
approximating  a  least  squares  fit  to  the  experimental  data  in  the  post- 
reflected  pressure  zone.  The  method  of  least  squares  is  a  numerical  analy¬ 
sis  technique  for  selecting  a  particular  curve  to  fit  some  given  data.  When 
approximating  a  least  squares  fit,  the  approach  is  to  maintain  an  area  under 
the  fitted  curve  that  equals  the  area  under  the  experimental  curve.  The 
applicability  of  this  approach  to  interpreting  shock  tube  test  data  is 
verified  in  Figures  11  and  12  in  Section  I I -2 ,  Error  Effects. 

Large  data  fluxuations  In  the  experimental  shock  wave  data  dictate  the 
need  for  an  approximated  least  squares  fit  Instead  of  a  computational  fit. 

A  major  characteristic  of  the  actual  computational  method  of  least  squares 
is  that  It  puts  great  emphasis  on  large  fluxuations  and  little  emphasis  on 
small  fluxuations.  As  a  result,  extreme  fluxuations  in  the  recording  of 
data  usually  dominate  the  results. 

Figure  3  is  a  typical  plot  of  experimental  post-reflected  data 
fitted  with  the  correspondi ng  approximated  least  squares  curve.  Note 


the  fitted  post-reflected  curve  follows  the  general  path  of  the 
experimental  data,  while  ignoring  large  fluxuations.  Ideal  post-reflected 
curves  characteristically  exhibit  a  steady  decline  in  pressure  with  time; 
however,  experimental  curves  may  decline  more  erratically  as  a  result  of 
test  conditions.  Whatever  the  post-reflected  curve  profile  may  be,  post- 
reflected  peaks  (Pr)  always  occur  at  the  initial  stagnation  time  (ts)  of  the 
post-reflected  zone. 

Figures  4  through  6  represent  panels  that  have  ideal  stagnation  times 
(ts)  of  approximately  .004  seconds  (according  to  Glasstone  calculations). 
Interpretations  of  Figures  4  through  6  yield  experimental  stagnation  times 
(ts)  between  .004  seconds  and  .007  seconds.  Ideal  stagnation  times  were 
used  as  guidelines  to  predict  ranges  where  experimental  stagnation  times 
should  occur. 

reflected  post- reflected 
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Figure  3  -  Approximated  least  squares  fit  on  a  typical  plot  of  experimental 
post-reflected  data. 


II-l.C  SPIKE  PEAKS 


Experimental  spike  peaks  are  dependent  upon  their  associated  post- 
reflected  pressure  peaks.  Relationships  between  spike  peaks  and  unreflected 
peaks  are  developed  in  detail  in  reference  2.  This  section  develops  the 
general  applications  of  these  ideal  relationships  as  applied  to  experimental 
pressure  data. 

Reflected  pressure  spikes  are  characteri stic  of  shock  waves  traveling 
non-parallel  paths  relative  to  the  surface  which  they  strike.  Ideal  spike 
peaks  for  reflected  shock  waves  that  strike  at  normal  incidence  to  a  flat 
surface  are  given  by: 


7P0  ♦  P 

Where:  Pr  =  Reflected  spike  peak  (psi) 

P0  =  Ambient  pressure,  ahead  of  the  shock  front  (psi) 

P  =  Peak  incident  overpressure,  behind  the  shock  front 
(psi). 

Peak  Incident  overpressures  were  read  from  the  experimental  data 
labeled  Tunnel  Wall  Incident  Overpressure  Time  History.  It  was  found  that 
substituting  post-relected  peaks  (Ps)  for  peak  overpressures  (P)  yields 
accurate  results  for  reading  the  experimental  data. 


Table  1  lists  Pr  to  Ps  relationships  within  the  range  of  the  experimen¬ 
tal  data.  These  relationships  are  used  to  approximate  spike  peak  magnitudes. 


Specific  magnitudes  are  determined  by  the  data  profiles  within  the  approxi 
mated  regions.  Examples  of  spike  peak  readings  are  given  in  Figures  4 
through  6. 


ps  Cpsi) 

Pr  (psi) 

.i 

2.01  x  Ps 

i 

2.06  x  Ps 

5 

2.28  x  Ps 

10 

2.53  x  Ps 

Table  1  -  Pr  to  Ps  relationships  for 
approximating  spike  peak  magnitudes. 


II-l.D  SPIKE  PEAK  WIDTHS 


Reflected  pressure  spike  peak  widths  are  determined  by  experimental 
data  profiles  at  the  spike  peak.  Ideal  blast  waves  do  not  exhibit  spike 
peak  widths;  however,  test  conditions  can  Induce  this  phenomenon.  Spike 
peak  magnitude  and  width  are  the  most  important  data  profiles  to  be  read, 
since  they  initiate  the  greatest  structural  and  material  responses.  Figures 
4  through  6  exhibit  data  taken  from  various  experimental  plots. 

Corresponding  notations  define  the  approach  applied  in  interpreting  both 
spike  peak  magnitude  and  width. 


II-l.E  INTERPRETATIONS 
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Figure  4  -  An  example  of  data  interpretation. 


Curve  Fitting  for  Figure  4 

Apply  a  least  squares  approximation  to  curve  fit  post-reflected  data. 
Note  that  experimental  data  fluxuatlons  subside  around  .028  seconds,  giving 
an  Indication  of  where  the  optimum  approximated  pressure  levels  are  for  the 
fitted  curve.  The  post-reflected  peat  (Ps)  Is  Interpreted  to  be  0.5  psl  at 
an  experimental  stagnation  time  (ts)  just  over  .006  seconds. 

Table  1  Indicates  the  reflected  peak  (PP)  Is  approximately  1.2  psl. 
Therefore,  the  Initial  experimental  peak  Is  taken  as  Pr  at  1.3  psl  and  the 
second  experimental  peak  Is  disregarded  since  It  Is  well  above  the  expected 
range  of  Pr.  Fitting  only  one  point  In  the  range  of  Pr  results  In  the 
absence  of  a  spike  peak  width. 
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Figure  5  -  An  example  of  data  interpretation. 


Curve  Fitting  for  Figure  5 

Approximate  a  least  squares  curve  to  fit  post-reflected  data. 
Experimental  data  fluxuatlons  subside  around  0.022  seconds  and  0.03  seconds, 
giving  an  Indication  of  optimum  fitted  curve  pressure  levels.  There  are 
several  data  fluxuatlons  before  0.016  seconds  that  are  Ig.iored.  The  post- 
refkcted  peak  (?,_.  is  Interpreted  to  be  approximately  .65  psl, 
corresponding  to  an  experimental  stagnation  time  (ts)  around  .005  seconds. 

Calculating  the  reflected  peak  yields  Pf  approximately  equal  to  1.3  psi 
-  refer  to  Table  1.  The  Initial  experimental  peak  Is  recorded  as  Pr  at  1.4 
psl  and  the  second  experimental  peak  Is  disregarded  since  It  Is  not  within 
the  expected  range  of  PP.  There  Is  no  spike  width  because  only  one  point  is 
fitted  In  the  range  of  PP. 
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Figure  6  -  An  example  of  data  interpretation. 


Curve  Fitting  for  Figure  6 

Curve  fit  the  post-reflected  data  with  a  least  squares  approximation. 
Post-reflected  data  clearly  deviates  from  theory  between  ts  and  0.012 
seconds.  This  is  shock  tube  phenomenon  and  is  recorded  as  fitted  data  since 
project  objectives  are  to  validate  NASTRAN  against  shock  tube  test.  Near 
0.032  seconds  experimental  data  fluxuatlons  subside,  giving  an  indication  of 
optimum  approximate  pressure  levels  for  the  fitted  curve.  The  post- 
reflected  peak  (Ps)  is  interpreted  to  be  approximately  2.3  psi, 
correspondi r.g  to  an  experimental  stagnation  time  (ts)  .0045  seconds.  Note 
that  Ps  is  defined  as  the  initial  pressure  of  the  post  reflected  curve. 

Table  1  approximates  the  reflected  peak  Pr  at  4.5  psi.  Experimental 
data  contains  two  points  at  this  pressure  range;  therefore,  a  spike  width 
does  exist  as  Indicated  by  the  fitted  curve.  This  deviates  from  the  ideal 
but  is  recorded  in  order  to  duplicate  shock  tube  phenomenon.  Smooth  spike 
peak  widths  are  fitted  and  experimental  data  fluxuatlons  at  the  peak  are 
disregarded. 


1 1-2  ERROR  EFFECTS 


Analytical  deflections  are  dependent  upon  interpretations  of  the 
experimental  reflected  pressure  time  histories.  The  i nterpretation  methodo¬ 
logy  developed  in  Section  I I -1  Is  subject  to  variations  from  one  interpreter 
to  another.  Considering  these  variations,  the  following  study  was  made  to 
gain  some  insight  to  the  degree  of  error  Induced.  This  study  consists  of 
two  approaches  as  follows:  spike  peak  width  variations,  and  complete 
displacement  of  the  pressure  time  histories.  Results  of  this  study  are 
discussed  In  Section  IV. 

Spike  peak  width  variations  were  analytically  applied  to  a 
22''x22"x.l93"  flat  unstiffened  panel  with  2  sides  clamped  and  2  sides 
pinned.  A  pressure  time  history  was  developed  for  the  first  shot  and  the 
spike  peak  width  was  altered  for  the  subsequent  shots  two  and  three.  These 
curves  were  developed  solely  for  the  purpose  of  observation  and  do  not 
necessarily  represent  Ideal  pressure  time  histories  as  developed  by 
reference  2.  The  specific  pressure  time  histories  developed  are  listed  with 
their  corresponding  plots  In  Figure  7.  Resulting  deflections  are  plotted  in 
Figure  8. 

Complete  displacement  of  pressure  time  histories  were  studied  to 
observe  the  effects  of  general  variations  In  the  interpretations  of  Iden¬ 
tical  shock  tube  blasts.  Actual  interpretations  were  developed  by  the 
methodology  covered  in  Section  1 1 -1 .  Ideal  interpretations  were  developed 
from  experimental  data  measured  in  the  regions  of  the  shock  tube  walls. 

These  curves  were  developed  for  the  sole  purpose  of  providing  various 


Interpretations  for  identical  shock  tube  blasts,  and  do  not  necessarily 
represent  ideal  pressure  time  histories  as  developed  by  reference  2.  Two 
comparisons  were  made  on  22"x22"x.l92"  flat  unstiffened  panels,  one  with  all 
sides  clamped  and  another  with  all  sides  hinged.  The  ideal  verses  actual 
interpretations  are  listed  with  their  correspondi ng  plots  in  Figures  9  and 
11.  Resulting  deflections  are  plotted  in  Figures  10  and  12. 


Figure  8  -  Error  Effects  Study  #1.  Deflections  resulting  from  pressures 
depicted  in  Figure  7. 


Figure  9  -  Error  Effects  Study  #2.  Complete  displacement  interpretation 
variation. 


IDEAL  vt.  ACTUAL  DEFLECTIONS 


Ill  NASTRAN  ANALYSIS 


A  total  of  seven  NASTRAN  models  were  developed  to  simulate  structural 
response  to  shock  tube  overpressures.  NASTRAN  models  were  developed  and 
compared  against  shock  tube  tests  performed  upon  four  flat  panels,  one  flat 
stiffened  panel,  and  one  honeycomb  panel.  One  curved  panel  was  studied 
qualitatively  since  insufficient  test  data  was  provided  for  a  quantitative 
study.  The  four  flat  panels  and  the  curved  panel  were  constructed  with 
CQUAD2  elements.  For  the  flat  stiffened  panel,  CQUAD2  elements  were  used 
for  the  skin  and  CBAR  elements  were  used  for  the  stiffeners.  CQUAD1  ele¬ 
ments  were  used  to  construct  the  honeycomb  panel .  Refer  to  "The  NASTRAN 
User's  Manual"  for  detailed  explai nations  of  these  elements.  The  seven 
NASTRAN  models  developed  are  presented  in  this  section.  Refer  to  Table  2 
for  general  model  specifications.  Figures  13  through  15  for  model  geometry. 
Criterion  for  building  NASTRAN  models  is  simplicity  of  design.  This  assures 
that  NASTRAN' s  effectiveness  will  be  researched  from  both  aspects  of  economy 
and  accuracy. 

Panel  deflections  are  used  as  the  criteria  for  comparisons  between 
NASTRAN  analysis  and  shock  tube  data.  Stress  was  not  used  as  a  criterion 
because  deflection  data  in  reference  one  is  of  more  consistent  quality. 
Deflection  and  stress  exhibit  a  linear  relationship  in  a  material's  elastic 
range.  Therefore,  deflection  is  a  valid  measure  of  NASTRAN's  effectiveness 
for  predicting  sure  safe  panel  response.  This  is  discussed  in  further 
detc.'i  in  Sectic  IV,  Discussion.  Deflections  are  coiroared  at  panel 
centers.  Tables  3  through  9  list  the  interpreted  pressure  time  histories 
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for  each  panel  analyzed.  Figures  18  through  24  plot  the  corresponding 
deflections  for  each  panel.  Results  of  the  comparison  are  discussed  in 
Section  IV. 
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Figure  13  -  NASTRAN  finite  element  model  for  Panels  1  through  5.  Refer  to  Table  2  for  specific 
panel  characteristics. 
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Figure  15  -  NASTRAN  finite  element  model  for  Panel  7,  a  curved  panel.  Refer  to  Table  2  for  specific 
panel  characteristics. 
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PANEL  6 

STIFFENED  CROSS  SECTION 


NASTRAN  versus  shock  tube  deflections. 
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NASTRAN  versus  shock  tube  deflections 


Figure  23  -  NASTRAN  versus  shock  tube  deflections. 
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IV  DISCUSSION 


Results  of  the  NASTRAN  analysis  and  shock  tube  test  comparisons  are 
listed  in  Table  10.  Results  of  the  error  effects  studies  are  listed  in 
Table  11  through  13.  First  deflections  characteristically  exhibit  the 
largest  deflection  responses  for  aperiodic  loading;  therefore,  magnitudes 
of  first  deflections  are  the  comparison  criteria.  Times  of  deflections  do 
not  dictate  stress  levels  and  are  therefore  considered  insignificant. 

Percent  error  between  magnitudes  of  first  deflections  was  the  measure  of 
effectiveness  in  both  the  NASTRAN  versus  shock  tube  comparison  and  the  error 
effects  studies. 


A  hand  calculation  shows  the  correlation  between  deflection  and  stress 
response.  The  calculation  determines  edge  stress  for  Panel  1  from  NASTRAN' s 
predicted  deflection  and  compares  this  to  edge  stress  data  measured  during 
the  shock  tube  test.  Equations  are  taken  from  Formulas  For  Stress  and 
Strain,  5th  Ed.,  Raymond  J.  Roark  and  Warren  C.  Young. 


For  rectangular  plates,  all  edges  fixed, 

uniform  load  over  entire  plate:  ^EDGE  =  <1  b2 

t£ 


Where:  (3^  =  0.3078 


a  =  0.0138 


ZMAX  =  cv  q  b4 


Constants  for  a  square  plate: 
Aspect  ratio  =  1 .0 


*•  •/  v  V  V 
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N-vy; 


q 


Uniform  static  load 


b  =  22.0  in  -  long  edge,  all  edges  the  same  for  panel  1 

t  =  .192  -  plate  thickness 

E  =  11.0  x  10^  -  Youngs  modulus 

0EDGE  =  Maximum  stress  at  edge 
^MAX  =  *219  in  (first  deflection  maximum). 

First,  calculate  the  equivalent  static  to  dynamic  uniform  load: 

9  =  ZMAX  E  t3 

O'  b* 

q  =  (.219  in)  (11.0  x  106  psi)  (.192  in)3 

(0.0138)  (22.0  in)* 


q  =  5.27  psi 

Second,  calculate  maximum  stress  at  the  edge: 

O’EDGE  ~  q  -- 
t2 


Oedge 


.3078)  (5.27  psi)  (22.0  in)2 
(.192  in)2 


^EDGE 


21297.18  psi 


Shock  tube  test  eoge  stress  measuremer.t  for  panel  1  show  a  maximum  value 


aEDGE 


25000  ps-; 


%  ERROR  =  NAS TRAN  -  Test  x  100 
I  Test  I 

%  ERROR 

%  ERROR  *  14.8 

This  corresponds  to  a  12.4%  error  in  the  deflection  comparison  for  panel  1. 
Therefore,  stress  analysis  does  correlate  very  closely  with  deflection  ana¬ 
lysis,  as  is  expected  since  stress  and  deflection  exhibit  a  linear  rela¬ 
tionship  within  a  materials  elastic  range. 

Stress  analysis  with  NASTRAN  requires  the  appropriate  model.  Such  a 
model  should  incorporate  a  center  element  and  refined  elements  at  the  middle 
of  the  longest  fixed  edge  of  the  panel.  A  center  element  is  required  to 
calculate  stress  at  the  panel  center.  Refined  elements  at  the  middle  of  the 
longest  fixed  edge  are  required  to  calculate  the  maximum  stress  for  a  fixed 
edge  panel.  The  refinement  of  elements  is  necessitated  by  the  sharp  stress 
gradient  that  occurs  at  a  panel's  fixed  edge.  Deflection  models  do  not 
require  such  element  refinement,  and  therefore  require  less  computer  time 

than  stress  models.  For  these  reasons  deflection  was  used  as  the  criterion 

* 

for  comparing  NASTRAN  results  with  the  shock  tube  tests  measurements. 

Inherent  errors  of  interpretations  of  pres,  re  time  histories  taken 
from  Boeing  Military  Airplane  Company  shock  tube  test  data  are  a  source  of 


21297.18  -  25000 
25000 


x  100 


error  in  the  NASTRAN  analysis  comparison.  Section  1 1-2  studies  two  possible 
error  effects.  First,  a  study  was  conducted  to  observe  the  error  effects  of 
spike  peak  width  variations.  Second,  a  study  was  conducted  to  observe  the 
error  effects  of  completely  displacing  the  pressure  time  history. 

Table  11  list  first  and  second  deflections  for  shots  1,2,  and  3  of  the 
spike  peak  width  study.  The  term  shot  refers  to  a  pressure  time  history. 
Figure  7  plots  the  three  shots  and  list  their  correspond!- ng  pressure  time 
histories.  Corresponding  deflection  data  is  plotted  in  Figure  8.  These  are 
considered  reasonable  variations  of  interpretations  for  spike  peaks  repre¬ 
sented  by  the  Boeing  Wichita  shock  tube  test  data. 

Six  possible  error  effects  are  taken  from  this  study  and  the  results 
are  listed  in  Table  12.  The  procedure  of  this  study  observes  each  shot  as 
an  actual  and  measures  the  error  effect  of  the  two  subsequent  shots  as 
ideals.  Results  of  this  study  show  that  it  is  reasonable  to  expect  approxi¬ 
mately  20%  error  from  a  spike  peak  and  37.4%  error  in  a  worse  case.  It  is 
emphasized  here  that  not  all  of  THUNDERPIPE1  s  pressure  time  history  data  is 
subject  to  sucn  inherent  interpretational  error. 

Complete  displacement  of  the  pressure  time  history  curve  is  the  second 
error  effect  study.  Figures  9  through  12  plot  the  pressure  time  histories 
av  resultant  deflections  of  the  two  cases.  Table  13  list  magnitudes  for 
the  first  and  second  deflection  peaks  and  their  relative  percent  errors. 
Times  of  deflections  are  not  listed  sine-  comparisons  are  made  against  iden- 
t'.c*  JAStkA*  ..10  i,  which  results  in  identical  times  v  reflections. 


Each  case  represents  two  interpretations  for  an  identical  shock  tube 
test.  Methods  of  interpretations  are  covered  in  Section  III.  Observing 
first  deflections,  a  39.2%  error  is  found  in  the  worse  case.  While  this 


does  not  represent  a  reasonable  inherent  interpretational  error  of  the 
Boeing  Military  Airplane  Company  shock  tube  data,  it  is  noteworthy  in  that 
it  emphasizes  the  effect  of  the  spike  peak.  Specifically,  it  takes  relati¬ 
vely  large  variations  in  interpretations  of  complete  displacements  to  pro¬ 
duce  the  equivalent  error  resulting  from  small  variations  in  interpretations 
of  spike  peak  characteristics. 

Results  of  the  comparison  between  the  NASTRAN  analysis  and  shock  tube 
test  are  listed  in  Table  10.  Figures  18  through  24  plot  the  correspondi ng 
deflections.  As  aforementioned,  percent  errors  between  magnitudes  of  first 
deflections  are  the  measurement  criteria  of  effectiveness. 

Panels  1  through  4  exhibit  very  close  comparisons  between  NASTRAN  data 
and  shock  tube  test  data  —  ranging  from  6.5  to  15.3%  error.  These  four 
panels  are  flat  and  homogeneous.  They  differed  in  aluminum  alloy,  panel 
thickness,  boundary  constraints,  and  pressure  time  histories.  Shock  tube 
test  data  for  these  four  panels  were  well  defined  by  Boeing  Military 
Airplane  Company,  and  therefore  considered  to  be  correctly  modeled  by 
NASTRAM . 

Trends  for  error  due  to  modeling  techniques  cannot  be  deduced  by  com¬ 
paring  these  four  panels.  By  relating  panel  descriptions  in  Table  2  to 
relative  pc -cent  errors  in  Table  13,  it  is  deter::.  ,  .vsd  that  neither  boundary 
constraints,  panel  thickness,  or  material  propertied  are  proportional  to 
magnitude  of  error. 
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Panel s  5  and  6  are  flat  nonhomoyeneous  panels.  Panel  5  is  a  honeyconb 
construction  and  panel  6  is  a  stiffened  panel .  Shock  tube  test  data  for 
panel  5  was  well  defined  by  Boeing  Wichita,  and  therefore  considered  to  be 
correctly  modeled  by  NASTRAN.  Accordingly,  panel  5  exhibits  a  very  close 
comparison  at  7.2%  error.  Panel  6  exhibits  the  worst  case  for  deflection 
comparisons  at  42.69%  error. 

Factors  that  may  have  affected  the  results  of  panel  6  are:  incorrect 
boundary  conditions,  inherent  interpretational  error  of  the  pressure  time 
history,  and  exceeding  the  linear  analysis  capabilities  of  NASTRAN .  This 
panel  was  modeled  with  boundary  constraints  as  stated  in  the  Boeing  Wichita 
final  report.  Since  deflection  frequencies  between  the  NASTRAN  and  shock 
tube  test  data  coincide,  it  is  assumed  that  boundary  conditions  are  defined 
reasonably  well.  Inherent  interpretational  errors  of  pressure  time 
histories  have  been  addressed  in  the  study  on  error  effects  and  show  that 
considerable  error  can  be  induced.  NASTRAN  uses  linear  finite  element  ana¬ 
lysis,  making  it  reliable  in  the  elastic  range  of  a  material's  response. 
Table  8  shows  that  panel  6  was  subjected  to  a  maximum  reflected  overpressure 
of  4.75  psi.  Plastically  yielding  deflections  during  shock  tube  test  exhi¬ 
bit  substantially  larger  deflection  magnitudes  than  NASTRAN,  since  NASTRAN 
continues  linear  past  the  yield  point  on  the  stress  strain  curve.  Boeing 
Military  Airplane  Company  documents  that  panel  6  plastically  deformed  during 
four  shock  tube  tests.  Figure  23  indicates  that  plastic  deformation  may 
have  occured  during  this  shot;  verifying  the  possibility  that  NASTRAN' s 
elas.ic  limits  mo/  nave  been  exceeding. 


Panel  7  is  used  for  qualitative  comparisons  only,  since  there  is  no 
experimental  deflection  data  available.  It  is  a  curved  homogeneous  panel 
which  was  subjected  to  a  maximum  overpressure  of  10.5  psi  --  the  largest  of 
all  panels  studied.  Figure  24  shows  the  deflection  response  predicted  by 
NASTRAN.  The  magnitude  of  the  first  deflection  is  relatively  small,  at  .043 
inches,  for  the  size  of  reflected  pressure  experienced.  This  coincides 
reasonably  with  Boeing  Military  Airplane  Company  documentation  that  panel 
seven  exhibited  no  permanent  deformation  after  the  test. 
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1st  and  2nd  Peak  Deflections 


Panel/ 

Defl. 

Mag 

nitude  (inches) 

Time  (sec) 

NASTRAN 

Test 

%  Error  (1) 

NASTRAN 

Test 

%  Error  (1) 

1  /1st 

.219 

.250 

12.4 

.0035 

— 

.0030 

16.7 

l/2nd 

-.088 

-.160 

45.0 

.0070 

.0065 

7.7 

2/lst 

.173 

.150 

15.3 

" 

.0060 

.0044 

36.4 

2/2nd 

-.040 

-.065 

38.5 

.0130 

.0090 

44.4 

3/lst 

.113 

.130 

13.1 

.0040 

.0040 

0.0 

3/2nd 

-.059 

-.067 

11.9 

.0085 

.0080 

6.3 

4/lst 

.029 

.031 

6.5 

.0025 

.0030 

16.7 

4/2nd 

-.012 

-.009 

33.3 

.0050 

.0055 

9.1 

5/lst 

.154 

.166 

7.2 

.0018 

.0032 

43.8 

5/2nd 

-.039 

-.060 

35.0 

.0036 

.0057 

36.8 

(  6/lst 

.066 

.115 

42.6 

.0025 

.0032 

21.9 

6/2nd 

-.042 

-.030 

_  .  .. 

40.0 

.0060 

.0055 

9.1 

7/lst 

.043 

(2) 

— 

.0025 

(2) 

— 

7 /2nd 

l - 

-.015 

.0045 

— 

(1)  %  Error  = 


NASTRAN  -  Test 
Te's't 


x  100 


(2)  Test  Data  Not  Available 


Table  10 
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Spike  Peak  Spike  Peak 

Width  Effects  Width  Effects 

(1st  and  2nd  Deflection  Peaks)  (1st  Deflection) 


Shot/ 
Def  1 . 

Time 

(sec) 

Defl . 

(i nches 

Shot 

Ideal /Actual 

%  Error  (1) 

1/lst 

.0039 

.1412 

1/2 

20.6 

l/2nd 

.0084 

-  .0572 

1/3 

27.2 

2/lst 

.0039 

.1779 

2/1 

26.0 

2/2nd 

.0084 

-.0928 

2/3 

8.3 

3/lst 

.0045 

.1940 

3/1 

37.4 

3/2nd 

.0086 

— 

-.1091 

3/2 

9.1 

Table  11  Table  12 


Complete  Displacement  Effects 
(1st  and  2nd  Deflections) 


Case/ 

Defl. 

Ideal 

Defl. 

(inches) 

Actual 

Defl. 

(inches) 

%  Error  (1) 

1/lst 

.2162 

.219 

1.3 

l/2nd 

-.0656 

-.088 

25.5 

2/lst 

.2408 

.173 

39.2 

2/2nd 

-.0739 

-.040 

84.8 

Table  13 


(1)  %  Error  =  Ideal  -  Actue '  I  x  ioo 
!  Actual 


V  CONCLUSIONS 


(1)  NASTRAN  is  an  accurate  analysis  code  for  predicting  elastic 
structural  response  to  shock  tube  tests  used  to  simulate  nuclear  blast 
overpressure  effects. 

(2)  Accurate  pressure  time  histories  of  shock  blast  are  extremely 
vital  for  accurate  predictions  of  structural  response. 

(3)  NASTRAN's  modeling  flexibilities  allow  for  greater  analysis  capa¬ 
bilities  of  nuclear  blast  overpressure  effects  than  are  allowed  with  present 
nuclear  effects  analysis  codes. 

(4)  NASTRAN's  programming  efficiency  results  in  less  computer  time 
required  than  with  present  nuclear  effects  analysis  codes. 

(5)  NASTRAN's  accuracy  in  overpressure  analysis  requires  accurate 
model  generation,  which  is  dependent  upon  accurate  structural  and  load  input 
data. 
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